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Abstract 

We compare four different methods for the numerical computation of the pure gluonic am- 
plitudes in the Born approximation. We are in particular interested in the efficiency of the 
various methods as the number n of the external particles increases. In addition we inves- 
tigate the numerical accuracy in critical phase space regions. The methods considered are 
based on (i) Berends-Giele recurrence relations, (ii) scalar diagrams, (iii) MHV vertices and 
(iv) BCF recursion relations. 



1 Introduction 



The fast and accurate computation of multi-parton amplitudes in QCD is essential for our under- 
standing of multi-jet processes at the LHC. It is a well-known fact that the conventional approach 
- summing up all Feynman diagrams - already reaches its limit for Born amplitudes when the 
number of external partons is about 5. More efficient methods construct the amplitudes recur- 
sively from smaller building blocks. The recursive approach allows to recycle information from 
already calculated pieces. An example of these methods are the Berends-Giele recurrence rela- 
tions [1-5]. The basic building blocks here are currents with one parton off-shell. 

Recently, a number of new methods for the calculation of helicity amplitudes in QCD have 
been introduced, motivated by the relationship of QCD amplitudes to twistor string theory [6]. 
In the Cachazo - Svrcek - Witten (CSW) construction [7, 8], tree level QCD amplitudes are 
constructed from vertices that are off-shell continuations of maximal-helicity-violating (MHV) 
amplitudes [9], connected by scalar propagators. Subsequently a set of recursion relations has 
been found [10-13] that involve only on-shell amplitudes with shifted, complex external mo- 
menta. In [14, 15] a method has been presented, which is close in spirit to the Berends-Giele 
recursion relations, but which involves only a set of primitive vertices and scalar propagators. 

These new methods enriched our understanding of the structure of QCD amplitudes. In 
particular they gave a precise answer to the question of what complexity we should expect in the 
result for a particular helicity amplitude, if we go beyond the simple maximal-helicity-violating 
ones. As an application, the singular behaviour of tree amplitudes in the multi-collinear limit 
could be derived from these methods [16, 17]. Furthermore these new methods turn out to be 
very useful in the construction of one-loop amplitudes [18-32] in conjunction with the unitary- 
based method [33,34]. 

Besides this undeniable progress in the understanding of the analytical structure of QCD am- 
plitudes, one question immediately arises: Do the new methods lead to improved algorithms for 
the computation of Born amplitudes in a purely numerical approach ? To answer this question we 
examine in this paper four different methods for the numerical computation of the pure gluonic 
amplitudes in the Born approximation. The methods considered are based on (i) Berends-Giele 
recurrence relations, (ii) scalar diagrams, (iii) MHV vertices and (iv) BCF recursion relations. 
We compare the efficiency of the various methods as the number n of external particles increases. 
In addition we investigate the numerical accuracy in critical phase space regions. 

This article is organised as follows: In section 2 we present the four different methods for 
the computation of the pure gluonic amplitudes in the Born approximation. In section 3 we 
compare the performance of the various methods and study the numerical stability in critical 
phase space regions. Section 4 contains our conclusions. Technical details are collected in the 
appendix. Appendices A and B define our conventions for spinors and Feynman rules. Appendix 
C collects the g — > gg splitting functions, which describe the collinear limit. In appendix D we 
comment on various optimisation techniques, which can be used to speed up the computation. 
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2 Description of the different methods 



The tree level amplitude with n external gluons may be written in the form 

*„(*?*,...,**») = g n - 2 £ 2Tr(r-(i)..T^))A„(^,...,^), (1) 

ceS„/z„ 

where the sum runs over all non-cyclic permutations of the external gluon legs. The symbol kj 
denotes the four-momentum of the j-th gluon and Xj its helicity. g denotes the strong coupling 
constant and T a the colour matrices, which are normalised such that Tv{T a T b ) = 1/2 h ab . The 
quantities A n , called the partial amplitudes, contain the kinematic information. They are colour- 
ordered, i.e. only diagrams with a particular cyclic ordering of the gluons contribute. 

In the computation of observables and cross-sections, the amplitude squared enters, summed 
over all helicities and colour degrees of freedom. We have 

\M 2 = 2>-y-%; £ £ |A„(^»...,^')f + °(^). CO 

Xu...,X„aeS n /Z n VW 

where N c = 3 denotes the number of colours. The colour-suppressed terms consist of interference 
terms between partial amplitudes with different colour-orderings. 

In this paper we investigate various methods for the efficient computation of the partial am- 
plitudes A n . We would like to separate this issue from the sum over all colour structures implicit 
in eq. (2). Therefore we focus on the quantity 



Xi,...,X„ 



(3) 



M n gives the leading-colour contribution to eq. (2), but in the context here it should be regarded 
as a quantity which depends only on the kinematical information and which helps to study the 
efficiency of various methods to calculate the partial amplitudes A n . It should be noted that only 
half of the helicity configurations need to be calculated, since parity relates a partial amplitude 
to the one with all helicities reversed: 

A n (k-i 1 , k n n ^j = —A n {k^ l ,...,k n . (4) 

In addition we investigate the numerical accuracy of the various methods in critical phase space 
regions. These are regions where one or more partons become unresolved. The simplest case 
involves single unresolved configurations, where one parton becomes either soft or collinear to a 
second parton. In the limit where one gluon j becomes soft, the partial amplitudes behave as 

A B+ i(*i,...*+ ..,*„ + i) v^^Jfc^|_A B (*i,. .,*„+!), 

— \ kj soft j— [kj+lkj—l] 

A n+ i(ki,...kj ,...,k n+ i) — > [k - + ik ■] [k k •••^ w + 1 )- ( 5 ) 
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(kikj) and [kikj] denote spinor products, which are defined in appendix A. The quantity M n 
factorises in the soft limit as 

M n+l (k h ...kj,...,k n+1 ) ^ 2 {2kj - 1 ^ - M n (k u ...,k n+l ). (6) 

(2kj-ikj){2kjk j+ i) 

In the collinear limit tree-level partial amplitudes factorise according to 

I Split_,(^,^)A„(...,^,...), (7) 

x=± 

where k a and k\, are the momenta of two adjacent legs. In the collinear limit we have K = k a + kb, 
k a = zK and kb = (1 — z)K. X, X a and Xb denote the corresponding helicities. The splitting 
functions are listed in appendix C. In the collinear limit the quantity M n behaves as [35-37] 

k a \\h 2 ( 2 , 2 A \ _ , ^ , , 8 



where the spin-correlation is given by 

S n = £ |£A B (^ 1 ,..,^,..,^-)+£*A B (*}' 1 ,..,^-,..,^-) 



(9) 



and 



£ = z <^ + l^+>. (10) 

2.1 Berends-Giele type recurrence relations 

Berends-Giele type recurrence relations [1,2] build partial amplitudes from smaller building 
blocks, usually called colour-ordered off-shell currents. Off-shell currents are objects with n 
on-shell legs and one additional leg off-shell. Momentum conservation is satisfied. It should be 
noted that off-shell currents are not gauge-invariant objects. Recurrence relations relate off-shell 
currents with n legs to off-shell currents with fewer legs. 

The recursion starts with n=\: 

J*(ki) = eP(k h q). (11) 

£^ is the polarisation vector of the gluon and q an arbitrary light-like reference momentum. We 
have the explicit formulae 

mK lH > V2{q-\k+y " K V2(k+\q-) 
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The recursive relation states that a gluon couples to other gluons only via the three- or four-gluon 
vertices : 



<M l^n\ — 

K 2 



1,71 



7=1 



n—2 n—l 



j=U=j+l 



,(13) 



where 

«ij = *,• + *,•+! + ... + *,• (14) 
and V3 and V4 are the colour-ordered three-gluon and four-gluon vertices 

The gluon current J M is conserved: 



yf po = /(2^ W -^V°-^V P ). (15) 



XX U, = 0- (16) 



v !=l 

Therefore terms proportional to ■ and proportional to Kj +l n can be dropped in eq. (13) and, 
using momentum conservation, it is legitimate to use the slightly simpler expression 

V^{hMM) = i(g vp (k 2 -k 3 y + 2 g P%-2g^k p 2 ). (17) 

for the three gluon vertex in eq. (13). 

The partial amplitude A n {k\ l , ■ ■■,k„ n ) is obtained from the gluonic current J M (k\ x , ...,k^"~i) by 
multiplying by the inverse gluon propagator and contracting with the polarisation vector for 
gluon n: 

A n (k\\...^) = %{k n ,q)-(iKl^ l )ji>$\...fc{). (18) 

A close inspection of the recursion relation eq. ( 1 3) shows that only the quantities ' , . . . , kj j ~l ) 
which respect the original order need to be calculated. Therefore an efficient implementation 
stores a list of four-momenta 

[k h k 2 ,...,k n ] (19) 

and a list of helicities 

(20) 

in memory and passes to the subroutine just two integers i and j, indicating that the quantity 

j^,...,k%-;) (2D 

should be computed. 
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2.2 Recursive calculation with scalar diagrams 

A modification of the Berends-Giele recursion relation was advocated in refs. [2] and [14]. In 
this approach all summations over Lorentz indices are replaced by a sum over the two physical 
polarisations. This reduces the number of multiplications needed for a contraction from four 
to two. The resulting recurrence relation consists of scalar propagators and a set of primitive 
vertices. 

Let q be a null-vector, which will be kept fixed throughout the discussion. Using q, any 
massive vector k can be written as a sum of two null- vectors k ? and q [38]: 



(22) 



Obviously, if k 2 = 0, we have k = k* . Note further that 2kq = 2$q. Using eq. (22) we may 
associate a massless four-vector $ to any four-vector k. Using the projection onto $ we define 
the off- shell continuation of Weyl spinors as 



\k±) - \k b ±). 
(k±\ -> (k b ±\ 

We are going to use the following abbreviations: 



(23) 



ki \kj-\- 



m 

(i-\j±k\l- 



+ \kj 



H)±H 



k\- 



(24) 



In spinor products, the projections $ are always used. Let us define an "off-shell amplitude" 

O n (^k^jk^ 2 , ...,k^j , (25) 

depending on n external momenta ki and helicities A, ( -. The momenta need not be on-shell, but 
momentum conservation is imposed: 



E*7 = o. 

7=1 



(26) 



By definition, the off- shell amplitudes O n are calculated from all Feynman diagrams contributing 
to the cyclic-ordered partial amplitude A n , by using the off-shell continuation eq. (23) for all 
external spinors and polarisation vectors, and by using the axial gauge for all internal gluon 
propagators. Compared to off-shell currents, which are used in recurrence relations of Berends- 
Giele type, an off-shell amplitude may have more than one leg off-shell. By construction, if all 
external particles are on-shell, the off-shell amplitude O n coincides with the physical amplitude 
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A n . We have the following recurrence relation: 

B (#,...,#)= £ iVs^Vi,^!,^) (27) 

+ I E E ^fe-i.4-i.<i*)^- ^ *H'(-^-irt 

X,X',X"=±j=2l=j+l A 1J-1 

A ;,z-i A /,n-l 
where we define the two-point amplitude to be the inverse propagator: 

2 (kf,-Kfj) = -4 (28) 

The partial amplitude A n coincides with O n , if all gluons are on-shell: 

A n (k^, ...,k n n ) = O n {k^ 1 , ...,k n n ). (29) 

There is only a limited number of non-zero vertices, which are listed in appendix B.2. This 
allows for a high degree of optimisation in the calculation of these vertices. The double and 
triple sums over the intermediate helicities in eq. (27) reduce in all cases to three non-vanishing 
terms. 

On the other hand it should be pointed out that in this approach the four-valent vertices 
depend (as do the three-valent vertices) on the momenta attached to these vertices. This should 
be compared to the standard Feynman rules, which enter the Berends-Giele recurrence relation, 
where the four-gluon vertex in eq. (15) is independent of the momenta. 

As in the Berends-Giele recurrence relation, an efficient implementation stores the sequence 
of four-momenta and helicities in a central place and just passes two integers i and j to the 
implementation of the recurrence relation, indicating the starting and ending points. 



2.3 Recursive calculation with MHV vertices 

In the Cachazo - Svrcek - Witten (CSW) construction [7], tree level QCD amplitudes are con- 
structed from vertices that are off-shell continuations of maximal-helicity-violating (MHV) am- 
plitudes, connected by scalar propagators. In maximal-helicity-violating amplitudes all gluons 
except two have the same helicity. Compact formulae for these amplitudes have been known 
for a long time [9]. Using the off-shell continuation eq. (23) the MHV-amplitudes serve as new 
vertices: 

(JkY 



v(k + k- k- k+) - ifviY' 2 {jk) 

V n {K l , ...,Kj , ...,K k , ...,K n ) — lyVZJ 
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Each MHV vertex has exactly two lines carrying negative helicity and at least one line carrying 
positive helicity. 

Bena, Bern and Kosower [39] derived a recursive formulation, which allows to obtain vertices 
with more gluons of negative helicity from simpler building blocks: 

V n (k x l \...,jfr) = —^—£ E ^V {M+2)modn (k%...rf l ,(-K hl )-) 

\nneg ^) j =U=j+1 J^jj 

X V (j-l) mod n (tf+l , • • • , fc/Ti , ( -K(l+ 1) ; ( j- 1) ) + ) , (31) 

where n neg is the number of negative helicity gluons. The recursion starts if n neg is less than two. 

For n neg = or n neg = 1 the quantity V n (k l \ . . . ,k„") vanishes. For n neg = 2 it is given by eq. 
(30). Again, the partial amplitude A n coincides with V n , if all gluons are on-shell: 

A n {k\\...^) = V n {k\\. (32) 

There are two points which should be noted: First of all, there is a double sum in eq. (31), which 
over-counts each contribution (n neg — 2) times. This over-counting is compensated by the explicit 
factor 1 / (n neg — 2) in front. 

Secondly, it is no longer possible to work with a static list of four-vectors and helicities, as was 
the case for the first two methods. The recursion relation eq. (31) inserts the four-momenta —Kj t 
and —Kn+i\(j-i) into the cyclic order. Therefore the lists of momenta and helicities have to be 
updated at each step of the recursion. This is best implemented by a double-linked list, which 
allows for the insertion of the new elements without copying the remaining ones. 

2.4 Recursive calculation with shifted momenta 

Britto, Cachazo and Feng [10] gave a recursion relation for the calculation of the n-gluon ampli- 
tude which involves only on-shell amplitudes. To describe this method it is best not to view the 
partial amplitude A n as a function of the four-momenta kj, but to replace each four- vector by a 
pair of two-component Weyl spinors. 

In detail this is done as follows: Each four- vector K M has a bispinor representation, given by 

K AB = K V°AB> K V = 2 Ka ^ A - ^ 

For null- vectors this bispinor representation factorises into a dyad of Weyl spinors: 

k^ = & k A s = k A k^. (34) 

The equations (33) and (34) allow us to convert any light-like four-vector into a dyad of Weyl 
spinors and vice versa. Therefore the partial amplitude A n , being originally a function of the 
momenta kj and helicities Xj, can equally be viewed as a function of the Weyl spinors k J A , k^ and 
the helicities Xf. 

A n (k^\ ...,kn") = A n (k A ,k^,'ki, ...,k A ,k^,'k n ). (35) 
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Note that for an arbitrary pair of Weyl spinors, the corresponding four- vector will in general be 
complex- valued. If (A, n ,A,i) 7^ (+, — ) we have the following recurrence relation: 



A n (k A ,k l B ,Xi,...,k A ,k n B ,X n ) = (36) 

n—l, . . 

E E A j (kA^ k l)M,kl,klM, ■■■,k J A ^\k J r\lj_ l JK A JK 6 ,-U 

j=3%=± 

x — 2 A M _y+2 (Ka,K b , X, k^k^Xj, ...,k A ,k B ,'k n -i,k A ,k^,'k n S j . 



If (X„,Xi) = (+, — ) we can always cyclic permute the arguments, such that (X M ,Xi) 7^ (+, — ). 
This is possible, since on-shell amplitudes, where all gluons have the same helicity, vanish. In 
eq. (36) the shifted spinors k l A , k n B , Ka and K B are given by 



kA = k A -zk\, K A = 



^(l + \K\n+y 



where 



k b- V {l + lK]n+ y (37) 



^45 = E 44' 4,-1 = det (38) 



and 



(39) 



(l + \K\n+) 

The recurrence relation starts with n = 3. The only non- vanishing amplitudes are 

A 3 {k A ,k B ,-,kl,kl,-,k A ,k B ,+) = i ^Jy^{73)^\) ' 

A 3 (*i,*J,+,44»+»44»-) = fv ^ [32] [ pi][13] ' (40) 

plus the ones with cyclic permutations of the helicities. It should be noted that due to the partic- 
ular choice of shifting the spinors, the three-point function with k A vanishes if the helicities are a 
cyclic permutation of (— , — , +). Similarly, the three-point function involving k n B vanishes if the 
helicities are a cyclic permutation of (+, +, — ). To speed up the computation the Parke-Taylor 
formulae in eq. (30) may be used for n>4. 

As in the previous method we have to update at each step in the recursion the list of Weyl 
spinors and the helicities. 
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n 


4 


5 


6 


7 


8 


9 


10 


11 


12 


Berends-Giele 


0.00005 


0.00023 


0.0009 


0.003 


0.011 


0.030 


0.09 


0.27 


0.7 


Scalar 


0.00008 


0.00046 


0.0018 


0.006 


0.019 


0.057 


0.16 


0.4 


1 


MHV 


0.00001 


0.00040 


0.0042 


0.033 


0.24 


1.77 


13 


81 




BCF 


0.00001 


0.00007 


0.0003 


0.001 


0.006 


0.037 


0.19 


0.97 


5.5 



Table 1 : CPU time in seconds for the computation of the n gluon amplitude on a standard PC (2 
GHz Pentium IV), summed over all helicities. 



3 Performance and numerical stability 

3.1 Performance 

We have implemented all four methods into numerical programs. For an unbiased comparison of 
the efficiencies of the different methods, each author has programmed all four methods indepen- 
dently, in order to eliminate possible dependencies on the programming skills of the programmer. 
It turned out that all programs gave the same pattern in the study of efficiency and accuracy. 

All methods give identical results within an accuracy of 10~ 12 for randomly chosen non- 
exceptional phase space points and up to 12 external particles. To investigate the performance in 
terms of CPU time we study the quantity <M n defined in eq. (3): 



M n = £ |a b (**V.,*J") 

X\,...,X n 



(41) 



It is clear from the algorithms that the first two methods (Berends-Giele and scalar diagrams) 
need a constant amount of CPU time for each helicity configuration, whereas the last two meth- 
ods (MHV and BCF) are very efficient if the helicities are predominately all plus or all minus, 
but take more CPU time if the helicity configuration contains roughly the same number of plus 
and minus helicities. To compare the different methods, the quantity M n sums over all helicity 
configurations. This corresponds to the situation encountered in the calculation of cross-sections 
and observables. Table 1 shows the CPU time needed for the computation of M n as n varies from 
4 to 12. The test was done on a standard PC with a 2 GHz Pentium IV processor. 

As can be seen from the table, the Berends-Giele type recurrence relation is the fastest 
method, as the number of external gluons increases. In second place comes the method with 
scalar diagrams. As already discussed in the presentation of the algorithms, these two methods 
are fast due to the fact that they can work with a static list of four- momenta and helicities. This 
avoids copying large amounts of data at each step of the recursion. The scalar diagram technique 
allows for a higher degree of optimisation in the subroutines, but this is out-weighted by the fact 
that in the Berends-Giele method each three- or four-valent vertex is called exactly once, whereas 
in the scalar diagram method each vertex is called three times with different helicity configura- 
tions. Table 2 shows the timings for the Berends-Giele method and the scalar diagram method 
for the computation of M n as n varies from 13 to 20. It should be noted that for n = 20 the results 
of the two methods agree within 10~ n . It can be seen from tables 1 and 2 that the time required 
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n 


13 


14 


15 


16 


17 


18 


19 


20 


Berends-Giele 


2 


4 


11 


27 


64 


149 


367 


831 


Scalar 


3 


6 


15 


36 


85 


195 


465 


1043 



Table 2: Continuation of table 1 for n in the range from 13 to 20 for the Berends-Giele method 
and the scalar diagram method. The settings are as in table 1 . 



by the scalar diagram method grows slower than that required by the Berends-Giele method as 
the number n of external particles increases. 

The MHV method is rather slow compared to the other three methods. This is related to the 
double sum appearing in eq. (31), which explicitly over-counts each contribution. In addition, 
the look-up tables we used to speed up the calculation are in this case rather memory-intensive. 
That is the reason why we were not able to compute the 12-point amplitude within this approach. 

The BCF method is faster than the Berends-Giele method as long as the number of external 
particles is below 9. For applications to three- or four-jet rates at LHC the BCF recurrence 
relations are therefore an improvement in efficiency. 

3.2 Numerical stability 

We have already mentioned that all methods give identical results for randomly chosen non- 
exceptional phase space points within an accuracy of 10~ 12 . In this section we study the nu- 
merical stability near exceptional phase space points, e.g. near singular configurations where 
one or more partons become unresolved. We limit ourselves to single unresolved configurations, 
where an external momentum becomes soft, or two external momenta become collinear. In these 
cases the quantity M n exhibits an infrared singularity and factorises into a singular function and 
a lower-point amplitude, as described by the equations (6) and (8). The singular behaviour can 
cause problems with the numerical stability of amplitude calculations. To investigate this prob- 
lem, we evaluated M n for configurations approaching each kind of singular limit. To illustrate 

(f) 

the stability, we have plotted in Figures 1 and 2 the ratio of to its factorised form as 
given by the right-hand sides of eqs. (6) and (8). The soft limit (Fig. 1) is described by x — > 
where x is the fraction of the total energy carried by the soft gluon. The onset of instability 
is at x ~ 10 -12 . The collinear limit (Fig. 2) is described by pr/E — > where pj is the trans- 
verse momentum involved in the collinear splitting (with z= 1/2). Instability occurs when the 
dimensionless variable pr/E ~ 10~ 7 . We observe no significant differences between the four 
methods. 

In addition to these physical singularities, spurious singularities might occur. An example 
can be found in the on-shell recursion relation. The shift in the spinors introduces sandwiches 
of the form (p,- — |/?& + pi\pj—) in the denominator. For example, the analytical formula for the 
six-gluon partial amplitude A6(l + ,2 + ,3 + ,4~,5~,6~) reads: 

A 6 (l+,2+,3+4-,5-,6-) = 



11 



|M 7 /Mf » | 



A n 



■17.5 



■15 



]□□□□□□□□□□□□□□□□□□□□□□□□□□□ 1 



0.8 



■12.5 



■10 



-7.5 



. 6 



. 4 



0.2 



-2.5 



log 10 x 



Figure 1 : Ratio of the sum of squared helicity amplitudes to its factorised form for a set of 7- 
gluon configurations where one gluon becomes soft, x is the energy fraction of the soft gluon. 
Key: o Berends-Giele * scalar diagrams A MHV rules □ on-shell. 



Ai 



>-|l+2|3-> : - 



+ 



(4-|5 + 6|l-) ; 



(61) (12) [34] [45]s 126 (2 - 1 1 + 6|5— > (23) (34) [56] [61]si 56 (2 - 1 1 + 6|5 



. (42) 



This introduces unphysical singularities when sums of external momenta become collinear. Of 
course, these cancel exactly in the final result, but can lead to problems when the recursion 
relation is evaluated numerically. An example of this is shown in Figure 3. Here we consider an 
amplitude of the form shown in eq. (42), in the limit that p\ + p^ becomes collinear to P2 + Ps- 
We have plotted the fractional error in the on-shell results for the sum over helicity amplitudes 
by comparing to those of the Berends-Giele recursion relation. The onset of instability occurs 
when the transverse momentum is of the order of lO" 7 ^. 

The other recurrence relations can also exhibit spurious singularities, as each require an arbi- 
trary light-like "reference" vector q to be specified, and various quantities diverge if this vector 
becomes collinear to one of the external momenta. For the Berends-Giele recurrence relations 
this vector is needed to define the polarisation vectors in eq. (12), and for the scalar diagram and 
MHV vertex approaches it is needed to fix the on-shell projection in eq. (22). The dependence 
of our results on q as q becomes collinear to an external momentum k is illustrated in Fig. 4. 
The scalar and MHV results become unstable when \fLq ~ 10~ 7 E, whereas the Berends-Giele 
recurrence relation is stable down to y/Lq ~ 10 12 E. This behaviour is expected, since in the 
Berends-Giele recurrence relation the reference vector q enters only the external polarisation 
vectors, whereas in the other two methods it also affects the internal lines. 

Overall, all methods exhibit satisfactory numerical stability properties. As far as spurious 
singularities are concerned, the Berends-Giele method performs slightly better. 
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□□□□□□□□□□□□□□□□□ 



□□□□□□□□□□□□□□□□□□□ D 



|M 7 /Mf » | 
2 



1 . 75 
1 . 5 
1 .25 
1 

. 75 
. 5 
0.25 



-14 



■12 



■10 



log 10 (pi/E) 



Figure 2: Ratio of the sum of squared helicity amplitudes to its factorised form for a set of 7- 
gluon configurations where two gluons becomes collinear. pr/E is the transverse momentum of 
the pair of gluons, normalised to the total energy. Key: o Berends-Giele * scalar diagrams A 
MHV rules □ on-shell. 



fractional error 
. 1 



0.08 



0.06 



0.04 



0.02 



log 10 (px/E) 



Figure 3: Fractional error in the sum of squared helicity amplitudes computed with the on-shell 
recursion relations for a set of 6-gluon configurations where k\ + &2 becomes collinear to £3+^4. 
Pt/E is the transverse momentum between the 2 pairs of gluons, normalised to the total energy. 
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fractional error 
0.04 r 



. 03 







. 02 







.01 











: ' • • 1 ■ ..- .■ ■ — 1 — — — 1 — ' — ' — 1 — 1 — ■ — 1 — ' — 1 — - 
2 -10 -8 -6 -4 -2 



log 10 (Vk.q/E) 



-12 



-0.01 



Figure 4: Fractional error in the sum of squared helicity amplitudes as the reference vector (q) 
used in the definition of each recursion relation becomes collinear with an external momentum 
(k). Key: o Berends-Giele ★ scalar diagrams A MHV rules. 

4 Conclusions 

In this paper we studied numerical implementations of recursive methods for the computation of 
Born gluon amplitudes. These amplitudes (together with corresponding ones, where additional 
quarks or vector bosons are involved), are relevant for LHC physics. They enter numerical LO 
or NLO program codes. As these calculations are based on Monte-Carlo integration over the 
phase space, the efficiency of the computation has a direct impact on the running time of the 
Monte-Carlo program. 

From the four methods considered, we found the Berends-Giele method performs best, as 
the number of external partons increases (n > 9). However, for a not so large number of exter- 
nal partons (n < 9), the on-shell recursion relation offers the best performance. As this is the 
range most relevant to LHC physics, this new method leads to an improvement in the numerical 
computation of Born amplitudes. 

We also investigated the numerical stability and accuracy. Here, all methods give satisfactory 
results. 



A Spinors 



For the metric we use 



gpv = diag(+l, -1,-1,-1). 



(43) 
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We define the light-cone coordinates as 

P+=Po+P3, P-=Po~P3, P±=Pi+ip2, P_L*=Pi-ipi- (44) 

In terms of the light-cone components of a null-vector, the corresponding massless spinors (p ± | 
and \p±) can be chosen as 

e J 3 / - pi _ t \ e '3 / p+ 



P+) = ^ „ I \P~) 



p+ J V\pT\\pi 

e 2 e '2 

(P+l = -/j=f (-/>!»/>+). (*H = — 7j=f(p+>P±*)> (45) 

vIp+I V\p+\ 

where the phase (|) is given by 

P+ = \p+\<P- (46) 

Spinor products are denoted as 

(pq) = {p- \q+) = p A qA, [qp] = (q+\p~) = qxP A - (47) 
We will also use the notation 

(p±\k\q±) = (p±\k„f\q±). (48) 

B Feynman rules 

B.l Colour-ordered Feynman rules 

The Feynman rules for colour-ordered partial amplitudes read: 





A' 



(49) 
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B.2 Scalar diagrammatic rules 

The non- vanishing primitive vertices involving only gluons are: 



o 1 

o 
o 
o 



v 3 (k- 1 k 2l k+) 



iV2(12) 



[lq][2q] 



iV2 



(12) 4 
(12)(23)(31) : 




V 3 (kt,ki, k J) = «V5[21]^ = i^^P] ■ 




.1 + 



, [lg]ggK3g)(4g) / ( g _|2-3| g -)( g -|4-l| g -) ' 
i (l9)(2?)[3 g ][4 ? ] V (q-\2 + 3\q-)(q-\4 + l\q-). 



4- 1 + 

= Vz).(fcj ,^2 5^-3 5^4 ) 

, [lg] (2g) [3gj (4g) n q -\l-2\q-)(q-\3- 4| g-) 
<1«> [2«] (3*) [4*] V (9 - 1 1 + 2| <?-> <<? - |3 + 4| 



+<*- 


2-3|f-)(f- 


4-1 


9-) 




2 + 3|< ? -)(< ? - 


4+1 


9-) 




C Splitting functions 

In the collinear limit the all-gluon tree-level partial amplitudes factorise according to 

A n+l (...,k a ,k b ,-) k ^ £ SpKt_ x (&,kl b )A n (...,K\...), (51) 

X=± 

where k a and k b are the momenta of two adjacent legs. In the collinear limit we have K = k a + k b , 
k a = zK and k b = (1 — z)K. A,, X a and A,& denote the corresponding helicities. The splitting 
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functions are: 



spiit g+ (g + , g + ) = o, s P iit g -Gr,g-) = o, 

Split g+ (g+,g-) = v^-==^, Splits (*",*+) = -V2 



3 



3 




s P iv(g-,g + ) = - , s P iv(g + ,g-) = -V2- 



l 



(52) 



y/z(T=z)(ab)' 



D Optimising techniques 

In this appendix we comment briefly on optimisation techniques we employed to speed up the 
computation. These techniques apply to all the four methods discussed in the main text. Recur- 
sive formulations have the disadvantage that they evaluate the same quantity over and over again. 
To overcome this obstacle, look-up tables are employed. C++ already offers in the standard li- 
brary the data structure std: :map, which is well suited for look-up tables. The template takes 
the types of three classes as arguments 

std: :map< class_key, class_result, class_compare > look_up_table; 

where class_key represents a class, which stores the input information on which a result de- 
pends. class_result defines the data type of the results and class_compare is a class, which 
provides an operator 

class class_compare 
{ 

public : 

bool operator () (const class_key & a, const class_key & b) const; 



for comparing two instances of type class_key, which is needed to keep the look-up table 
sorted. The class std: :map needs of the order log(n) operations to look-up a specific entry, if 
the table is filled with n entries. This can be improved by using a method with constant look-up 
time. 

A second point concerns temporary variables. All methods operate on data structures, like 
four-vectors or spinors, which are rather time-consuming if they need to be copied. C++ allows 
operator overloading and for example the addition of fourvectors can be coded as follows 

fourvector p,q; 
fourvector sum = p + q; 

However, this first calculates the sum of p and q in a temporary variable and copies the result in 
a second step to sum. It is more efficient for selected operations to write a method 
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void add (fourvector & sum, const fourvector & p, const fourvector & q) ; 

which avoids copying temporaries: 

fourvector sum,p,q; 
add (sum, p, q) ; 

This is not as elegant as in the first code fragment, but more efficient when used in low-level 
routines. 
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